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ABSTRACT 

We examine the angular momentum loss and associated rotational spindown for mag- 
netic hot stars with a line-driven stellar wind and a rotation-aligned dipole mag- 
netic field. Our analysis here is based on our previous 2-D numerical MHD simula- 
tion study that examines the interplay among wind, field, and rotation as a function 
of two dimensionless parameters, one characterizing the wind magnetic confinement 
(r/» = Bl q Rl/ Mvoo), and the other the ratio (W = V rot /V or b) of stellar rotation to 
critical (orbital) speed. We compare and contrast the 2-D, time variable angular mo- 
mentum loss of this dipole model of a hot-star wind with the classical 1-D steady-state 
analysis by Weber and Davis (WD), who used an idealized monopole field to model 
the angular momentum loss in the solar wind. Despite the differences, we find that 
the total angular momentum loss J averaged over both solid angle and time follows 
closely the general WD scaling J = (2/3)MflR A , where M is the mass loss rate, Q. is 
the stellar angular velocity, and Ra is a characteristic Alfven radius. However, a key 
distinction here is that for a dipole field, this Alfven radius has a strong-field scaling 

Ra/R* ~ V* , instead of the scaling Ra/R* ~ ^/rj*~ f° r a monopole field. This leads to 
a slower stellar spindown time that in the dipole case scales as r spin = T mass l.bk/y/rj^, 

where T mass = M/M is the characteristic mass loss time, and k is the dimensionless 
factor for stellar moment of inertia. The full numerical scaling relation we cite gives 
typical spindown times of order 1 Myr for several known magnetic massive stars. 

Key words: MHD — Stars: winds — Stars: magnetic fields — Stars: early-type — 
Stars: rotation — Stars: mass loss 



1 INTRODUCTION 

In recent years, improvements in spectropolarimetry have 
made it possible to detect moderate to strong (10 2 — 10 4 G) 
magnetic fields in a growing number of hot, luminous, mas - 
sive stars of spectral type O and B (e.g.. lDonati et al] 2002). 
The high luminosity of such stars drives a strong stellar wind 
through line scatte ring of the star's continuum radiation 
l)Castor et alJ Il975l . hereafter C AK) . The first two papers 
in this series focus on developing numerical magnetohydro- 
dynamics (MHD) simulations of the confinement and chan- 
neling of this stellar wind outflow by a dipole magnetic field 
at the stellar surface. For m odels without rotation, Paper 1 
( ud-Doula & Owockilliool ) showed that the overall effect of 
the field depends on a single "wind magnetic confinement 
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parameter" rj, (defined in eqn. <(3j below) that characterizes 
the ratio of magn etic to wind energy de nsity near the stellar 
surface. Paper 2 ()ud-Doula et al.ll200ct l extended this study 
to include field-aligned rotation, examining a wide range of 
magnetic confinement parameters (rj* from near unity up to 
1000) in stars with equatorial speeds that are a substantial 
fraction (W =1/4, 1/2) of the critical (orbital) speed. Paper 
2 focussed mainly on formation and disruption of equato- 
rial rigid-body disks. The present study utilizes the same 
2-parameter MHD simulation study to examine the angu- 
lar momentum loss and associated stellar spindown for this 
case of wind outflow from a hot star with a rotation-aligned 
dipole. 

Most of the previous literature on magnetic wind spin- 
down has focussed on cool, solar-type stars, for which the 
wind is driven by the high gas pressure of a corona heated to 
more than a million K [see, e.g., reviews bv lMestell (|l968al lbl. 
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1 19841 ) and lMestel fc SpruitJ (|l987l l; iTout fc Pringld l| 19921 )]. 
The mechanical energy to heat the corona is thought to 
originate from the strong convection in the hydrogen re- 
combination layers below the stellar surface. Moreover, the 
interaction of this convection with stellar rotation is under- 
stood to drive a dynamo that generates a complex stellar 
magnetic field and activity cycle. Much of the emphasis in 
studying cool-star spindown has thus focussed on the feed- 
back of faster rotation in generating both a stronger field and 
more mechanical heating to drive a stronger coronal wind, 
which then act together to give a more abrupt spindown for 
younger, more-rapidly rotating stars. For a middle-age star 
like the sun, which is roughly halfway through its expected 
10 Byr main-sequence lifetime, the rotation speed is thus 
quite slow, only about 2 km/s at the solar equator, with an 
associated rotation period of about 27 days. 

For massive, hot, luminous stars with radiatively driven 
winds, the direct study and modeling of wind magnetic spin- 
down is more limited. This is partly due to a general expec- 
tation that the absence of a hydrogen recombination convec- 
tion zone means that hot stars should not have the magnetic 
dynamo activity cycles of cooler stars. Nonetheless, as noted 
above, spectropolarimetric observations have directly de- 
tected large-scale fields in a growing number of O-type (cur- 
rently 3) and early-B-type (ca. 2 dozen) stars, often charac- 
terized by a more or less constant dipole that is tilted relative 
to the stellar rotation axis. This steady nature and large- 
scale contrast with the active and complex magnetic activity 
cycle of cool stars, and suggests a primordial fossil origin in- 
stead of active dynamo generation. However there have been 
models based on active generation in the convective core 
llMacGregor fc Cassinelli 20031; MacGregor fc Charbonneaul 



ll999T l or envelope (|Mullan fc MacDonaldll200ll . I2005T I. 

All the directly detected, oblique dipoles with strong 
fields also exhibit a clear periodic rotational modulation in 
circumstellar signatures like X-ray emission and in wind sig- 
natures like the UV P-Cygni line profiles. But weaker, as-yet 
undetected, smaller-scale fields might also be one cause of 
less regular wind variations, such as the episodic discrete 
absorption components (DACs) commonly seen in the ab- 
sorption troughs of P-Cygni profiles. 

Compared to cooler stars, the rotation in early-type 
stars is quite fast with inferred periods typically one to sev- 
eral days, and projected surface rotation speeds (Vsini) of 
order one or two hundred km/s. This has even been used to 
argue for a lack of wind magnetic spindown, and thus for 
a lack of a d ynamically significant magne tic field. Based on 
estimates by iFriend fc MacGregorl (|l984 hereafter FM84) 
for the dependence of spindown ti me on hot-star m a ss loss 
rate and magnetic field strength, iMacGregor et all (|l992T l 
argued that magnetic fields in massive stars with still-rapid 
rotation must generally be less than about 100 G. 

In the context of the present paper, a key issue for this 
FM84 analysis of wind magnetic spindown in hot stars is 
that it is based on the idealization - fi rst introduced in the 
seminal paper bv lWeber fc Davisl (1 19671 . hereafter WD67) on 
spindown from the solar wind - that the radial field at the 
stellar surface can be described as a simple monopole. Al- 
though magnetic fields can never be actual monopoles, this 
idealization greatly simplifies the analysis, making tractable 
a quasi-analytic, 1-D formulation. Moreover, for the sun it 
is somewhat justified by the inference that, beyond a few 



solar radii from the solar surface, the outward expansion 
of the solar wind pulls the sun's complex multipole field 
into an open radial configuration, roughly characterized by 
a "split monopole", with opposite polarity on the northern 
vs. southern sides of a heliospheric current sheet. 

However, as noted above, the fields detected in hot stars 
are often inferred to be dominated by a large-scale, dipole 
component. For such cases, it is not clear how applicable the 
monopole-based analysis of WD67 and FM84 should be for 
estimating spindown rates. Indeed, while some magnetic O- 
stars have slow rotation perio ds (e.g. 15 days for 6 Ori C, 
and ca. 500 d for HD191612; iDonati et all (|2006r i), several 
early-B stars with very strong (many kG) magnetic fields 
still retain a quite rapid rotation (e.g . period of 1.2 days for 
a Ori E, see lGroote fc Hungedll982h . 

The central purpose of this paper is to use our previ- 
ous MHD simulation parameter study to examine the wind 
magnetic spindown of massive stars with a rotation-aligned 
dipole field. Compared to the 1-D semi-analytic studies for 
an idealized monopole field, the numerical simulations of 
dipole fields here requires a second spatial dimension (lati- 
tude), but the restriction to field-aligned rotation allows us 
to retain a 2-D axisymmetry in which the variations in az- 
imuth are ignorable. However, our '2.5-D' formulation still 
retains the crucial azimuthal components of the flow velocity 
and magnetic field. 

As discussed in detail in section 3 below, this param- 
eter study shows that the loss of angular momentum for 
this dipole case is a highly complex, time-variable process, 
characterized by a quasi-regular cycle of build-up and re- 
lease of angular momentum stored in the circumstellar field 
and gas. But quite remarkably, the time- and angle-averaged 
total angular momentum loss follows a simple scaling rule 
that is quite analogous to that derived by WD67, with how- 
ever a key dipole modification in the scaling of the associated 
Alfven radius relative to that applicable to the simple WD67 
monopole model. As discussed in section 4, this leads to a 
substantial general reduction in the spindown rate for this 
dipole case relative to that derived by FM84 for hot-stars 
using a monopole model. A concluding section 5 summarizes 
the results and outlines directions for future work. To lay the 
groundwork for interpreting the detailed numerical simula- 
tion results, the next section gives some essential background 
on the basic analytic scaling laws for angular momentum loss 
in the WD67 monopole model and some minor variants. 



2 BACKGROUND 

2.1 The Weber &: Davis monopole model for the 
solar wind 

An outflowing wind carries away angular momentum and 
thus spins down the stellar rotation. Winds with magnetic 
fields exert a braking torque that is significantly larger than 
for non-magnetic cases, due to the larger lever arm of mag- 
netic field lines that extend outward from the stellar surface. 

A seminal analysis of this process was carried out by 
WD67, who modelled the angular momentum loss of the so- 
lar wind for the idealized case of a simple monopole magnetic 
field from the solar surface. In terms of the surface angular 
velocity Q and wind mass loss rate M, a key result is that 
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the total angular momentum loss rate scales as 

j = ~ MUR\ , (1) 

with Ra the Alfven radius, defined by where the radial com- 
ponents of the field and flow have equal energy density. This 
can be intuitively interpreted as the angular momentum loss 
that would occur if the gas were kept in rigid-body rotation 
up to Ra, and then effectively released. But while helpful 
as a kind of mnemonic, it is not literally the case, since in 
fact, as WD67 emphasized (and is discussed further below), 
most of the angular momentum is actually lost via Poynting 
stresses of the magnetic field, and not by the gas itself. 

For any radius r, the energy density ratio between radial 
field and flow is given by 



rj(r) 



pv?/2 



(2) 



The latter equalities emphasize this energy ratio can also 
be cast as the inverse square of the Alfvenic Mach num- 
ber, Ma = v t /Va, where the radial Alfven speed, Va = 
Br/y/4-Kp, with p the wind mass density. The Alfven radius 
is then defined implicitly by ti(Ra) = 1. 

Using detailed flow solutions of the equations for a gas- 
pressure-driven solar wind, together with in situ measure- 
ments of the radial magnetic field near earth's orbit, WD67 
estimated the solar wind Alfven radius to be Ra ~ 24. 3Rq . 
With this extended moment arm, even the quite low so- 
lar wind mass loss rate implies a substantial spindown over 
the solar lifetime, providing a possible explanation of the 
slow solar rotation . Applications to other s olar-type stars 
(e.g. lMestellll968allbt iMestel fc Spruitlll987l ; lKaroIerlll988l ; 
iTout fc Pringiel 1992) have largely focussed on the potential 
feedback of rotation on the field strength and mass loss rate. 

But in the present context of hot-star winds, for which 
the mass loss rate is set near the surface by the physics of 
radiative driving, we can derive approximate explicit expres- 
sions in terms of fixed values for the equatorial field strength 
B eq at the surface radius R* , and for the wind mass loss rate 
M and terminal flow speed Voo ■ Specifically, following papers 
1 and 2, if we define here a wind magnetic confinement pa- 
rameter, 



V* = 



B eq R* 



then we can write the energy density ratio in the form 

- 2 i „, 



T)(r) = rj* 



r 
R * 



v T (r) (r/i?») 2 (l- R*/ r y 



(3) 



(4) 



where q is the power-law exponent for radial decline of the 
assumed magnetic field, and the latter equality assumes now 
a monopole field (q — 2), together with a canonical 'beta' 
velocity law with index (3, and terminal speed i>oo- 

For the monopole case and velocity index of either (3 = 1 
or f3 — 2, an explicit expression for the Alfven radius can be 
found from solution of a simple quadratic equation, yielding 



^ = 1/2 + Vv« + 1/4 ; /? = 



and 



V^h ; = 2. 



(•») 



(6) 



For weak confinement, -q, <C 1, we find Ra — + R*, while 
for strong confinement, 77* 3> 1, we obtain _Ra — * ^/rj7R,. 
Application in eqn. (JTJ) then gives an explicit expression for 
the angular momentum loss rate. 

Note that in the strong magnetic confinement limit 
w, 2> 1, the scaling Ra ~ y/rjZ implies that the angular 
momentum loss for this monopole model becomes indepen- 
dent of the mass loss rate, 
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7J» > 1 . 



(7) 



For a star of moment of inertia I = kMR% (with typically 
k ~ 0.1), the associated characteristic spindown time for 
stellar rotation is 



J 
j 



kMQRj 
lMQ,R\ 



2-k 
2Z 

V* 



^kMvoo 
Bi q Ri 



(8) 



where the third equality gives a scaling in terms of a charac- 
teristic mass loss time, r m(lss = M/M, and the last equality 
gives the mass-loss-independent scaling. Note that although 
we derived this scaling in the context of line-driven hot-star 
winds, we did not make any explicit assumptions about the 
wind-driving mechanism. As such, we can apply this even to 
gas-pressure driven winds. In particular, for the solar wind, 
in situ measurements near 1 au give a speed w 400 km/s 
and radial field strength B au w 5 x 10~~ 5 G, translating to a 
monopole field strength B eq w B a u(au/ Rq) 2 ~ 2.3 G at the 
solar surface. Using k = 0.059, eqn. (JSJ then gives a solar 
spindown time of ca. 8.6 Byr, comparable to the solar age 
of ca. 5 Byr. 



2.2 Angular momentum loss from gas vs. 
magnetic field 

Let us next consider general expressions for angular mo- 
mentum loss, comparing the contribution due to the gas vs. 
the magnetic field. For spherical coordinates representing 
radius r, co-latitude 0, and azimuth <f>, let subscripts denote 
associated components of the vector velocity (v r , ve, V4) or 
magnetic field, {B r , Bg, B^). We are interested in the an- 
gular momentum about the rotation axis. For the gas, the 
associated angular momentum per unit mass is given by 
the azimuthal speed times the distance to the rotation axis, 
v$r sin#. Multiplying this by the mass flux density pv r then 
gives the angular momentum flux within an element of area 
r 2 d(j>dp (defining p, = cosO). Upon integration over azimuth 
(assuming axisymmetry) , we obtain the latitudinal distri- 
bution of gas angular momentum loss at any radius and 
colatitude, 



d Jnas . sin 6 



(9) 



where m = 4ivpv r r 2 gives the local mass loss rate (which in 
general could vary in latitude, radius, or time). 

For the magnetic field, the angular momentum loss is 
proportional to the r, (f> component of the Maxwell stress 
tensor, 



B r B$ 
4tt 



(10) 



which represents the radial flux density of azimuthal mo- 
mentum. As before, multiplying by the axial distance r sinf? 
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converts this into an associated angular momentum flux, 
which upon azimuthal integration gives the latitudinal dis- 
tribution of magnetic angular momentum loss, 



dj n 



-r B r B<t,r- 



(11) 



2.3 J 



These expressions for loss of rotational angular momentum 
apply for any general magnetic field, including the rotation- 
aligned dipole model discussed in detail in §3. But to il- 
lustrate some characteristic properties, let us first examine 
them for the simple WD67 monopole field model. In this 
case, both and B^ scale in proportion to sinf, giving an 
overall latitudinal dependence with sin 2 8 — I — (i 2 . For a 
slow rotation case like the sun, B r , v r , p and rh are other- 
wise largely independent of latitude 0- As such, latitudinal 
integration (from \i — —1 to /i = +1) gives an overall an- 
gular momentum loss that is just a factor 2/3 smaller than 
computed from the WD67 equatorial analysis, 



2 • 2 



rB r Bsr\ 



(12) 



Since under the "frozen-flux" condition (applicable to 
ideal MHD) the local velocity vector is parallel to the local 
field, we can relate the equatorial B$ and v$ through, 



Br 



V4, 



(13) 



Defining then an equatorial specific angular momentum 
j eq = J e q/rh, we find that combining eqns. (|12[) and (|13p 
gives for the gas specific angular momentum, 



Jga 



rv<j> = 



3e 



. M% - Q.r 2 



(14) 



At the Alfven radius Ra, where the Alfvenic Mach number 
Ma = 1, the denominator vanishes; ensuring continuity thus 
requires that the numerator also must vanish at this point, 
which implies j eq = QRa 2 ■ This thus provides the basis for 
the key WD67 scaling cited in eqn. ([1]). 

The fraction of angular momentum carried by the gas 
at any radius is then given by 



Jeq 



rv<p 
VR A 2 



VrA/Vr 



- v tA Ra /v r r 2 



(15) 



where v r A = v r (RA)- In the spherical expansion of the 
WD67 model, a similar 2/3 latitudinal correction applies 
to both the gas and total angular momentum, and so eqn. 
(|15p also gives the spherically averaged gas angular momen- 
tum fraction, J ga s/ J = jgas/jeq- In particular, at large radii, 
note that this gas fraction of angular momentum becomes, 



J a 



J 



= 1 



VrA 



(16) 



The remaining fraction is carried by the magnetic field, 
Jmag/J = iva/«oo- Since typically v t a/Voo < 1, the WD67 



1 However, in monopoles winds with more rapid rotation, the 
field and outflo w both tend to become defl ected toward the ro- 
tation pole; see ISuess &: Nernevl [|l975h and IWashimi fc Shibatal 
jl993h . 



monopole field model thus predicts that most of the angu- 
lar momentum is lost via the magnetic field, not the gas. In 
their analysis of the solar wind, WD67 obtained an asymp- 
totic ratio of about 3:1 for angular momentum of field to 
gas. 

In the somewhat broader context of a monopole field 
in a wind with velocity parameterized by a standard 'beta' 
law, v(r)/voc 



1 - R,/rf, we find for /3 = 1, 



J a 



J 



R, 
Ra 



1 



(17) 



y/T). + 1/4 + 1/2 

Analogous, but more complicated expressions can be derived 
for other values of the index f3. At the Alfven radius, appli- 
cation of L'Hopital's rule in eqn. (|15|l gives for general (3, 



Jgas 



j tot \ A P + 2(Ra/R*-1)' 1 ' 

For example, note that for strong fields in the j3 = 1 case, 
the gas fraction of angular momentum at Ra is just half the 
asymptotic value. 



3 ANGULAR MOMENTUM LOSS FOR A 
ROTATION-ALIGNED DIPOLE 

3.1 MHD simulation parameter study 

While the above simple monopole model is convenient for 
analytic study, actual magnetic fields on the sun and other 
stars can be far more complex, often represented by many 
higher order multipoles. As a first step in extending the 
above spindown analysis to a more physically realistic mag- 
netic configuration, let us consider now the case of a dipole 
field with axis aligned with that of the stellar rotation. Rela- 
tive to a monopole, such an aligned dipole implies variations 
in a second spatial dimension, namely co-latitude 0, as well 
as radius r, but still retains the axisymmetry that allows 
neglect of variations in azimuth <j>. Nonetheless, as shown in 
papers 1 and 2, the competition between wind outflow and 
closed magnetic loops now leads generally to an inherently 
complex, time-dependent behavior that is not amenable to 
direct analytic study, but instead requires numerical simu- 
lation through the solution of the equations of magnetohy- 
drodynamics (MHD). 

The MHD simulations in papers 1 and 2 have exam- 
ined specifically the effect of dipole fields in hot, luminous, 
massive stars with radiatively driven stellar winds. Paper 
2 presented a detailed parameter study of the competi- 
tion among wind, field, and rotation as a function of two 
dimensionless parameters, namely the wind magnetic con- 
finement parameter 77* defined in eqn. @, and a rotation 
parameter W = V ro t/V or b, representing the ratio of the 
equatorial surface rotation speed to equatorial orbital speed 
V or b = v GM/R t . The analysis in paper 2 focussed par- 
ticularly on the accumulation of wind material into a dense 
equatorial disk, confined in nearly rigid rotation between the 
Kepler co-rotation radius Rk = W~ 2 ^ 3 R*, and the Alfven 
radius Ra- 

The remainder of this paper now uses this same param- 
eter study to analyze the angular momentum loss in this case 
of a radiation-driven wind from massive star with a rotation- 
aligned dipole field at the stellar surface. The reader is re- 
ferred to paper 2 for full details of the numerical method, 
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spatial grid, and assumed stellar and wind parameters. How- 
ever, for the convenience of the reader we briefly summarize 
these below. 

For all our calc ulations we use the ZEUS-3D 
|Stone fc Normanlll992l) numerical MHD code. Our imple- 
mentation here adopts spherical polar coordinates with 300 
grid points in radius r spaced logarithmically and 100 grid 
points in colatitude 6 with higher concentration of points 
near the magnetic equator. We also assume symmetry in 
the azimuth <j> direction. To maintain this 2.5D axisymme- 
try, we assume the stellar magnetic field to be a pure dipole 
with polar axis aligned with the rotation axis of the star. 

As in paper 1 and 2, we consider only the radial compo- 
nent of radiative force with assumed flow strictly isothermal. 
To avoid the effects of oblateness and gravity darkening, we 
limit ourselves to moderate rotation rates applied to a model 
with stellar parameters characteristic of £ Pup (see table 1 
in paper 1). 



3.2 J for dipole scaling of Alfven radius 

A key result of paper 2 (see eqn. 9 therein) was that in this 
case of a dipole field, the equatorial Alfven radius follows 
the approximate scaling, 

^ « 0.29 + (r h + 0.25) 1/4 , (19) 
it* 

which represents an approximate solution of the quartic 
equation that arises from requiring i](Ra) = 1 in a dipole 
(q — 3) model with a (3 = 1 velocity law [cf. eqn. (j4])]. 

Note in particular that, in the strong confinement limit 
77* 2> 1, the Alfven radius in this dipole case now has the 
scaling Ra/R* ~ r\V A , instead of the scaling, Ra/R* ~ ?7* 
of the monopole model. As we show below, this modified 
scaling of the Alfven radius in a dipole model has important 
implications for the associated scaling of angular momentum 
loss. 

To proceed, let us introduce a basic ansatz that the over- 
all angular momentum loss of this rotation-aligned dipole 
case can still be described in terms of the simple WD67 
expression of eqn. {l}, if one just uses the dipole-modified 
scaling (|19|l for the Alfven radius. Specifically, let us define 
a "dipole- WD" (for "dipole Weber-Davis" ) angular momen- 
tum loss as 

Jdwo = I MQR A = | MQRl [0.29 + (J?, + 0.25) 1/4 ] * . 

(20) 

Here the mass loss rate M and wind terminal speed Vao 
used to compute the magnetic confinement parameter 77, 
are those the star would have without a magnetic field, e.g., 
as set by the physics of radiative driving. Even without a 
field, there is, however, a modest dependence of the mass 
loss rate on rotation, found here from numerical simulations 
of non-magnetic cases (see figure 8 of paper 2) to give about 
a 10% increase in going from the W = 1/4 to the W — 1/2 
rotation case. With these mass loss rates, we use this simple 
analytic form (|20[) to scale the numerical MHD simulation 
results presented below. 



3.3 Standard model case: 77, = 100 and W = 1/2 

As in paper 2, let us focus first on a standard case with 
moderately strong magnetic confinement, 77* = 100, and 
with rotation at half the critical rate, W = 1/2. At any 
time snapshot of the time-dependent simulation, we can use 
eqns. @ and pip to compute, at each colatitude 6 and ra- 
dius r, the latitudinal distribution of angular momentum 
loss associated with the gas and field, with their sum thus 
giving the associated total loss, dj/d/j. 

3.3.1 Spatial distribution and time variation of dj/dfi 

For each of the same 6 time snapshots shown in figures 2 
and 3 of paper 2, figure [T] here presents plots of dj/dfi, with 
the colorbar normalized in units of the predicted dipole- WD 
scaling of eqn. I|20p . The changes among the panels empha- 
size the intrinsic time variability of the model, with intervals 
of nearly stationary confinement (upper row) punctuated by 
episodes of sudden magnetic breakout (lower row). 

Nonetheless, particularly during the confinement inter- 
vals, there is a clear characteristic pattern for the overall 
distribution of angular momentum loss. Near the surface, 
dj/dfi is essentially zero within the close-field equatorial 
loops, but this is compensated by a concentration of angu- 
lar momentum loss in the open-field regions at mid-latitudes 
that, as we show below, is contributed mainly by the mag- 
netic component. Moreover, further from the star, there is 
a broad latitudinal distribution from the magnetic compo- 
nent combined with an equatorial concentration from the 
gas component. 

The colorscale in figure illustrates this latitudinal dis- 
tribution and time variation of the gas, magnetic and total 
(gas + magnetic) angular momentum loss dj/d/i along both 
the inner and outer boundary for this standard model case. 
For the inner boundary, the broad white region at low lati- 
tudes emphasizes quite clearly now that the closed magnetic 
loops above the equatorial surface represent a kind of "dead 
zone" with little or no angular momentum loss. Instead, the 
mid-to-high latitudes of open field carry a strong concen- 
tration of the surface angular momentum loss. In contrast, 
at the outer boundary, there is a broad latitudinal distri- 
bution of angular momentum loss from the field, together 
with a equatorial concentration from the gas. Both the in- 
ner and outer boundary distribution of angular momentum 
also show a clear time variation associated with the ca. 1 
Msec cycle of confinement, buildup, and release of material 
trapped in closed magnetic loops. 

3.3.2 Time variability of latitudinally integrated J 

The line plots in figure [3] compare the time dependence of 
the latitudinally integrated angular momentum loss, J, for 
the gas, field, and total components, again evaluated at the 
inner (red curves) and the outer boundaries (blue curves). 
At the inner boundary, the gas component is negligible, with 
most of the angular momentum associated with the field. At 
the outer boundary, the gas component is highly variable, 
but still generally small compared to the field. The total 
angular momentum loss shows a nearly periodic variation, 
characterized by a gradual ramping up that ends in a sudden 
drop, with however a slight relative lag in the outer vs. inner 
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Figure 1. For the standard case of rj, = 100 and W = 0.5 at the same time snapshots as in figures 2 and 3 of paper 2, plots of the 
spatial variation of dj/dfi, with the colorbar normalized in units of the predicted dipole-WD scaling of eqn. J20P with an estimated 
Alfven radius Ra = 3.45-R, . 



variation. This lag reflects a cycle of storage and release of 
angular momentum within both the circumstellar field and 
gas. 

The colorscale plots in figure [4] show also the full radius 
and time variation for the gas, magnetic, and total angu- 
lar momentum loss, with the colorbar again normalized in 
units of the predicted dipole-WD scaling of eqn. (|20[> . The 
results show quite vividly the intrinsic time variability, par- 
ticularly for the gas component, which varies from intervals 
of little or no angular momentum loss, to a series of radially 



ejected streams, punctuated by a strong interval of loss dur- 
ing the magnetic breakout. The generally pale color reflects 
the fact that the overall level of gas angular momentum loss 
is a small fraction of the total expected from the dipole-WD 
scaling, particularly near the stellar surface. By comparison, 
the magnetic component is stronger and less variable. Over- 
all, we see that the total angular momentum loss varies by 
about 50% above and below the predicted dipole-WD value 
from eqn. ([20]) . 
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Figure 2. Plot of the latitude and time dependence of the gas and field components of dj/dfj, through the inner boundary (r = R„; 
bottom panels) and the outer boundary (r = 15R*; top panels), for the standard model case. The colorbar is again normalized by the 
dipole-WD scaling of eqn. H20I I, 
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Figure 3. Time dependence of latitudinally integrated angular momentum loss J through both the inner and outer boundary, computed 
for the gas (left panel), field (middle) and the total for gas + field (right), in the standard model case. 



3.3.3 Spatial variation of time- averaged J for gas vs. 
magnetic component 

To show more clearly the overall spatial variation, figure [5] 
plots the radial dependence of the time-averaged angular 
momentum loss for both the gas (blue curve) and magnetic 
field (red curve), along with the total loss (black curve). 
Note that this time-averaged total loss is nearly constant in 



radius, at a value that is remarkably close - about 90% - to 
the simple analytic dipole-WD predicted scaling! 

Moreover, much as found in the WD67 monopole model, 
the angular momentum loss by the gas increases with radius, 
but still remains everywhere relatively small compared to 
the magnetic component. Indeed, the dashed curves compare 
the corresponding WD scalings given by eqn. (|15[) , using the 
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Figure 4. plots of the full radius and time variation of J, again computed for gas, magnetic, and total components in the standard 
model, with the colorbar normalized in units of the predicted dipole-WD scaling of eqn. 1201 . 
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Figure 5. Radial variation of gas, field, and total angular momen- 
tum loss, again scaled by the total loss in the dipole-WD model. 
The solid curves show the time-averages from the standard nu- 
merical simulation model, while the dashed curves compare the 
corresponding WD scalings implied by eqn. H15I I, using the Alfven 
radius Ra from the dipole eqn. I|19| l assuming a /3 = 1 velocity 
law and n* = 100. 



Alfven radius Ra from the dipole eqn. (|19J1 . and assuming 
a f3 — 1 velocity law and rj, = 100. 

The good general agreement shows that, despite the 
complex time variation of the dipole case, the overall, time- 
averaged scaling of angular momentum loss can be quite 
well modeled through the simple monopole scalings of WD, 
as long as one just accounts for the different scaling of the 
Alfven radius with the magnetic confinement parameter n* . 



3.4 Dependence on magnetic confinement and 
rotation parameters, r), and W 

To build on this success in using the dipole WD scaling to 
characterize angular momentum loss in this standard case 
of moderately strong confinement and rotation, let us now 
examine how well this simple scaling agrees with the numeri- 



cal simulation results for variations of magnetic confinement 
parameter and rotation parameter W . 

For the rotation cases W — 1/4 and W = 1/2, the lower 
and upper panels of figure [6] compare the variation of to- 
tal, time-averaged angular momentum loss rate vs. rj, (on a 
log-log scale) for both the numerical simulations (triangles) 
and analytic form p0[) (squares). The overall agreement is 
remarkably good for both rotation cases, confirming that, 
quite independent of the rotation parameter W, this very 
simplified form (|20[) provides a good description of the scal- 
ing of the average angular momentum loss in this case of 
aligned dipoles. 

Note that the ordinate axes in figure [6] are labeled in 
CGS units computed for the specific stellar model used in 
the simulations. But these specific values are essentially ar- 
bitrary. For any star of interest, the appropriate physical 
values can be readily derived from the dipole-WD scalings 
in eqn. (|20|l . Indeed, the plot can likewise be characterized 
as giving the inverse of the spin-down time, which in the 
dipole-WD model has the specific scaling 

Tmass ~ [0.29 + {f), + 0.25) 1 / 4 ] 2 ' 

Figure [7] shows the fraction of the gas and magnetic 
components of angular momentum loss at the outer bound- 
ary for each model, plotted vs. magnetic confinement pa- 
rameter r/* (again on a log scale), for both the W — 1/4 
(left) and W = 1/2 (middle) rotation cases. For comparison, 
the right panel plots the corresponding large-radius scalings 
for the pure monopole model with a /3 = 1 velocity law, 
as given by eqn. (|17|l . For low and moderate magnetic con- 
finement, 77* < 30, there is good general agreement between 
the analytic scalings and the numerical simulation results; 
but for stronger confinement, the numerical results show the 
gas fraction reaching a minimum and then increasing with 
increasing ry». 

The reasons for this increase are not entirely clear, but 
likely could be the result of the channeling and confinement 
of wind gas into the equatorial, nearly rigid-body disk dis- 
cussed in paper 2. The magnetic torquing that spins this 
material up into a rigid disk represents a transfer of an- 
gular momentum from field to gas that has no parallel in 
the pure-outflow, monopole models of WD. Once sufficient 
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Figure 6. Time-averaged angular momentum loss for all the 
models (triangles), plotted vs. n* (on a log-log scale), for both 
the W = 1/4 (left) and W = 1/2 rotation cases. The squares 
compare the scalings predicted by the dipole-WD approximation 
[1201 



material accumulates in this disk, the outward centrifugal 
force overwhelms the inward confinement of magnetic ten- 
sion, leading to a breakout of this material that now carries 
outward a strong gas component of angular momentum loss. 



4 DISCUSSION 

4.1 Role of mass-loss "Dead Zone" 

A principal result of the above parameter study is that the 
overall level of angular momentum loss from an early-type 
star with a rot at ion- aligned dipole can be well described by 
the simple dipole-modified WD scaling given in eqn. (|20p . 
In this formulation, the mass loss rate and wind terminal 
speed used to compute the magnetic confinement parameter 
77* and associated spindown are those the star would have 
without a magnetic field, as set by the physics of radiative 
driving. 

As shown in figure 8 of paper 2, the actual net mass loss 



rates in this parameter study of dipole winds show a signif- 
icant decline with increasing confinement parameter n*, fit 
roughly by the scaling relation [given in eqn. (24) of paper 2] , 

Mb 



Mb 



R,/R c + 1 - a/1 - 0.5 * R,/R K , (22) 



where R c w i?» + 0.7(i?A — R*) is a maximum "closure" 
radius of magnetic loops, and Rk = R*/W 2 ^ is the Ke- 
pler co-rotation radius. The former accounts for the effect 
of the mass loss "dead zone" of closed magnetic loops, while 
the latter corrects for the eventual centrifugal breakout that 
can occur from some initially closed loops above the Kepler 
radius. 

In previous discussions of rotational spindown of mag- 
netic winds, this dead zone has gener ally been presumed 
(e.g., lMestel|[l968al ; iDonati et all 12006) to lead to a down- 
ward modification in the net angular momentum loss that 
would otherwise occur, based on the notion that the mass 
trapped in these closed loops does not (at least for loops 
closing below the Kepler radius) escape from the star, and 
thus should not contribute to the angular momentum loss. 

This notion seems partly based on the perception that 
the gas itself is the principal direct carrier of the angular 
momentum loss. But both the WD67 analysis and the sim- 
ulations here show that the dominant effect of the gas is 
indirect, inducing an azimuthal field component that then 
carries the bulk of the angular momentum loss, particularly 
near the stellar surface. Figure shows that, even for this 
magnetic component, the closed loops at low latitudes do 
indeed represent a dead zone for loss of angular momentum, 
as well as mass. The net effect, however, seems merely to 
shunt a fixed total amount of angular momentum towards 
the mid-latitudes, carried by the azimuthal twisting of the 
open magnetic field. As the wind expansion opens up the 
field beyond the Alfven radius, this transport of angular 
momentum spreads to cover all latitudes (cf. right vs. left 
panels of figure[5J), and includes an increasing (but generally 
still minor) component for the gas (see figures [3] [4] [5] and 

So an important lesson of the above parameter study is 
that this additional dead-zone reduction in the net J does 
not apply to the dipole-modified WD scaling form (I20|l . In a 
sense, it is already incorporated in the reduction associated 
with the change from the monopole scaling J ~ Mn» to 
the dipole scaling J ~ M^/rj^. To see this, note that, if we 
ignore the minor rotational correction of the Kepler term, 
the net mass loss reduction given by eqn. (|22[) has the strong- 
confinement (??* ^> 1) scaling 



Mb 
Mb=o 



R, 
2R C 



R: t 



IARa 



Ra > R* 



(23) 



If we now use this to apply a "dipole dead-zone" correction 
to the standard monopole scaling for angular momentum 
loss, we then find 



J 



Mb Ra 



M b =oRa ~ M B =o0> , 



(24) 



which thus reproduces the dipole scaling using the non- 
magnetic value for the mass loss rate! 

The upshot then is that the dipole-modified scaling (|20[l 
using the non-magnetic mass loss rate effectively already 
accounts for the dead-zone reduction of the actual mass loss. 
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Figure 7. Fraction of the gas and magnetic components of angular momentum loss, plotted vs. magnetic confinement parameter n* (on 
a log scale), for both the W = 1/4 (left) and W = 1/2 (middle) rotation cases. The right panel shows analytic scalings for the pure 
monopole model with a j3 = 1 velocity law, as given by eqn. dlTI . 



Table 1. Estimated Spindown Time for Selected Known Magnetic Stars. 



Star a 


M/M e 


-R* /-R© 


P (days) 


k 


M (10~ 9 M Q /yr) 


Voo (1000 km/s) 


B p (kG) 




Tspin (Myr) 


9 1 Ori C 1 


40 


8 


15.4 


0.28 


400 


2.5 


1.1 


15.7 


8 


HD191612 2 


40 


18 


538 


0.17 


6100 


2.5 


1.6 


7.6 


0.4 


C Cas 3 


8 


5.9 


5.37 


0.1 


0.3 


0.8 


0.34 


3200 


65.2 


a Ori E 4 


8.9 


5.3 


1.2 


0.1 


2.4 


1.46 


9.6 


1.4 xlO 5 


1.4 


p Leo 5 


22 


35 


7-47 


0.12 


630 


1.1 


0.24 


20 


1.1 



a References: iDonati et al.l l|2002l l: foonati et all l|2006l l: ^Neiner et all l|2003l l and ISmith ". 
^Kholtvein et al.l ll2007f i 



Bohlenderl l|2007h : ferticka et"all l|2006h : 



4.2 Spindown time 

From the above, the overall stellar spindown time is pre- 
dicted to follow the scaling in eqn. (I21[l , In the strong- 
confinement limit, this gives 




B eq R* V M 



1.1 x 10 yr 



M/R, 



B p /kG Mq/Rq V M. 



(25) 



Comparison with the monopole scaling JH} shows that the 
spindown is no longer independent of the mass loss rate, 
but now varies with its inverse square root. In addition, the 
dependencies on surface field and radius is now inverse linear 
instead of inverse square. 

In practice, application of this scaling relation requires 
observational and/or theoretical inference of the various pa- 
rameters. The last equality in (|25|l gives characteristic B2- 
star scalings for the mass loss, M_ 9 = Af/(lO _9 M /yr), 
and wind speed, vs = v oo /(10 8 cm/s). Note that the mag- 
netic field is now quoted as a surface value at the pole, B v , 
and is scaled in kG. This is a typical value for known mag- 
netic massive stars, as inferred from Stokes V measurement 
of photospheric lines with circular pola rization by the Zee- 
man effect (see, e.g.. lDonati et al]|2002l ). 

The mass to radius ratio M/R* can best be estimated 
from atmosphere models for the given spectral type, but 
generally for massive main sequence stars this should be just 



somewhat above the solar ratio. The moment of inertia con- 
stant k can be estimated from stellar structure models, and 
should be typically k ~ 0.1, perhaps somewhat higher near 
the zero-age main sequ ence (ZAMS ), and then decreasing 
by up to 50% with age (|Claretll2004l) Pl 

Perhaps the most difficult parameters to infer are those 
for the stellar wind. Fortunately, these enter only in propor- 
tion to the square root of the ratio of wind speed to mass 
loss rate. But because the actual values for both of these are 
likely to be strongly affected by the magnetic field, it may 
generally be better not to infer them from direct observa- 
tions for the actual star in question, but rather to use the 
inferred spectral type to apply observational or theoretical 
scaling laws for the values in similar non-magnetic stars. 

Table 1 lists spindown times based on estimated param- 
eters for a sample of known magnetic hot stars. The first 
two known mag n etic O-stars, 0\ Ori C (|Donati et al.ll2002h 
and HD 191612 (|Donati et al.ll2006l ), are both slow rota- 
tors, with periods of respectively 15 d and 538 d. 6 1 Ori C 
is thought to be quite young, less than 0.2 Myr, and so 
still on the ZAMS, while HD 191612 is thought to be more 
evolved, with an ag e of 2-3 Myr. Using parameters quoted in 
IDonati et al.l [2006L we infer corresponding spindown times 
of respectively ~8 Myr and ~0.4 Myr. This implies that 
magnetic wind braking seems unlikely to explain the slow 



2 Note that this presumes solid-body rotation; if the stellar enve- 
lope should decouple from the core, the effective k could be much 
lower (to account for the lower envelope mass). But a strong in- 
ternal magnetic field should generally be effective in enforcing 
near-rigid rotation in the interior. 
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rotation of 8 1 Ori C, but it does seem potentially relevant 
for HD 191612. Alternatively, magnetic effects during star 
formation could have lead to an initial slow rotation. 

For He-strong stars, a key object is the B2V star 
a Ori E. This has an estimated polar field strength of 
BkG ~ 9.6, and with remaining parameters as in table 1, 
we estimate a spindown time of ~1.4 Myr. As a main se- 
quence B2 star, its age is likely comparable to this. The 
rotation period, 1.2 d, is still quite short, about twice that 
for critical rotation, implying only a moderate net spindown 
since formation. 



4.3 Extension to non-aligned dipoles and higher 
multipoles 

However note that, as is typical of magnetic massive stars, 
most of the above cases exhibit dipole fields with a signifi- 
cant tilt angle to the rotation axis. MHD simulation of such a 
tilted dipole requires accounting for 3-D, non-axisymmetric 
outflow as the polar field sweeps around in azimuth. This 
remains a challenge for future studies. Without such simu- 
lations, we can only offer some general speculations on how 
such a tilt might affect the spindown. Generally, it seems 
that angular momentum loss should be modestly enhanced, 
because of the factor two stronger polar field. 

But another factor might be the open nature of this po- 
lar field. One might expect this to lead to a larger magnetic 
moment arm, perhaps even following the stronger monopole 
scaling, J ~ 77*, rather than the dipole form J ~ y/rj7- But 
the analysis in the Appendix suggests that such magnetic 

polar-axis fields in an oblique rotation case should actually 

1/3 

follow a weaker spindown scaling, J ~ <qj . This implies 
that, much as in the aligned-dipole case, the overall, spher- 
ically averaged loss rate for an oblique-dipole wind should 
still be dominated by the regions surrounding closed loops, 
with a net scaling that thus is similar to the aligned case. 

In some magnetic hot stars, the inferred field is man- 
ifestly non-dipole. For example, the B2IVp He-strong star 
HD 37776 (V901 Ori) has been inferred to have quadrupole 
field that dominates any dipole component, with peak 
strength of ~10 kG (Thompson fc Landstreetl 1 19851 ). The 
luminosity class IV would normally imply an evolved, main- 
sequence star, but the still moderately short rotation period 
of 1.5387 d seems to suggest little spindown. If we assume 
a somewhat larger radius and higher mass loss rate than 
a Ori E, so that the non-magnetic parameters in (|25|l may 
be near unity, then applying the inferred field of 10 kG gives 
a spindown time of ~1.1 Myr. 

Moreover, although detailed predictions must await 3-D 
MHD simulations of such a quadrupole (or higher multipole) 
case, one might infer that the steeper radial decline of the 
quadrupole field should lead to a weaker spindown. For a 
general field scaling with r~ q , with q = 3 for a dipole and 
q — 4 for a quadrupole, the ratio of magnetic to wind en- 
ergy should decline as r 2 ~ 2q , implying an Alfven radius that 
scales as Ra ~ r]t^ 2q 2 \ or Ra ~ rjl^ 6 for a quadrupole. In 
the strong-confinement limit, the expected spindown scaling 
should thus become 



T~spin,quad 



V 



\k_ 

1/3 



(26) 



Compared to a dipole of the same surface strength, the spin- 



down time for a quadrupole would thus be about a factor 
r/l^ 6 longer. This might seem like a weak correction, but for 
He-strong stars of spectral type B2, the low mass loss means 
that a surface field of 10 kG gives a confinement parameter 
of 77, « 10 7 , thus implying a factor ~10 times longer spin- 
down for a quadrupole vs. a dipole of same surface strength. 
In this context, the expected spindown time for HD 37776 
becomes of order ~10 Myr. 

On the other hand, the presence of a more complex 
field might make it easier for the wind to break open the 
magnetic flux associated with large-scale dipole component, 
and so allow a more extended magnetic moment arm. In 
principal, this could even lead to a stronger spindown effect. 
iMikulasek fc et all (|2008T l have recently reported a direct 
measurement of rotational spindown in HD 37776. They in- 
fer a 17.7 ± 0.7 second increase in the 1.5387 day period 
over a span of 31 years, translating to a spindown time of 
just 0.23 Myr! This is substantially shorter than the above 
dipole estimate from wind torquing. Indeed, the luminosity 
class IV would normally imply an evolved, main-sequence 
star, but the still moderately rapid rotation together with 
the very rapid spindown seem to require a very young age. 

Overall, the study of angular momentum loss in these 
more non-aligned dipole or higher-order multipole cases 
must await further, 3-D MHD simulation studies. 



4.4 No sudden spindown during breakout events 

In a few magnetic stars, there appears to be evidence for 
sudden change in rotation period. For the rapid l y rot ating 
Ap star CU Virginis (HD 124224). iTrigilio et all (20081 ) cite 
radio observations indicating changes in rotational phase 
over 20 years, identified with two discrete period increases 
of 2.18 s in 1984 and about 1 s in 1999. In terms of the 
still-rapid rotation period 0.52 days, the associated average 
spindown timescale over this 15 year timespan is quite short, 
only about 300,000 years. These sudden period changes 
could be associated with changes in the star's internal struc- 
ture llStepierilll99Sl ); but ITrigilio et al.l i|2008h suggest they 
might also be the result of a sudden emptying of mass ac- 
cumulated in the star's magnetosphere. 

The simulations here do indeed show repeated episodes 
of substantial emptying, but careful examination indicates 
that these do not lead to sudden jumps in the angular mo- 
mentum of the star itself. Both figures [5] and [3] show, for 
example, that angular momentum loss through the stellar 
surface varies smoothly through the cycle of magnetosphere 
build-up and release, and if anything is actually less dur- 
ing the sudden breakouts that characterize magnetospheric 
emptying. In terms of stellar rotation, the outward transfer 
of angular momentum to the magnetosphere occurs gradu- 
ally, as the magnetosphere fills up, effectively increasing the 
moment of inertia of the star+magnetosphere system. As 
such, when the emptying does occur, it merely represents the 
final escape for angular momentum that had already been 
lost to the star. Moreover, most of the angular momentum 
is not even contained in the trapped gas, but rather in the 
stressed magnetic field, which figure 3b shows varies quite 
smoothly from the stellar surface. 

Thus, even apart from questions of the magnitude of 
angular momentum loss needed to explain the average spin- 
down inferred for CU Vir, magnetospheric emptying does 
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not seem well-suited to explaining the claimed suddenness 
of rotation changes in this star. 



4.5 Comparison with previous spindown analyses 

Past studies of wind magnetic spindown have primarily fo- 
cussed on cool, low-mass, solar-type stars. The main analysis 
aimed specifically at spindown for hotter, higher-mass star, 
the main analysis was by FM84, who derived 1-D, steady- 
state solutions for the equatorial flow of a CAK-type line- 
driven wind from a rotating hot star with a Weber-Davis 
style monopole magnetic field. To ensure their steady so- 
lutions passed smoothly through the various flow critical 
points (i.e., those associated with Alfven, slow- mode, and 
fast-mode MHD waves, as well as the usual CAK critical 
point), they had to use a quite elaborate numerical itera- 
tion scheme; as such, they do not quote any simple scaling 
forms for the Alfven radius and associated spindown rate. 

Instead, their Table 1 lists numerical results for a set 
of 16 models with various rotation rates and field strengths, 
assuming fixed stellar and wind parameters for a typical O- 
supergiant. We find here that these results can be generally 
well fit (within ca. 10-15%) by the simple general monopole 
scaling for Ra, eqn. ©. 

However, FM84 quote spindown times as low as 
60,000 years, i.e., for their case 4, with field strength of 
1600 G, and mass loss rate of 5.32 x lO _6 M0/yr. The asso- 
ciated confinement parameter is 77* « 80. Since the angular 
momentum loss rate for an aligned dipole is a factor y/rj7 
lower than for a monopole, we see that a dipole with a sim- 
ilar surface field strength (1600 G) as FM84 model 4 would 
have about a factor 9 longer spindown time, or now about 
0.5 Myr. In contrast, their weaker-field model, e.g. case 1 
with only 200 G, has a confinement parameter 77* as 1.7, 
implying only about a 30% dipole increase above the 2 Myr 
spindown time quoted for this case by the FM84 monopole 
analysis. Overall, the dipole modifications in spindown rate 
found here suggest t hat the upper lim i ts on surfaced field 
strength inferred by iMacGregor et all (|l992f ) for rotating 
hot-stars are likely to be too low. 

For cooler, solar-type stars with coronal-type pressure- 
driven winds, the literature on wind magnetic spindown 
is m ore exten s ive, and includes b o th semi-analytic stud- 
ies (jOkamotol Il974l : iMestell 1 19841 ; iTout fe Pringlel " 



and numerical iter a tions o r simulations, dSakurail 
Washimi fe Shibatal 1 19931 ; iKeppens fe GoedbloedT 



Matt fe Ba lick 2004). In addition to spindown by a mag- 
netized coronal wind during the star's main sequence phase, 
there has been extensive study of angular momentum loss 
during pre-main-sequence accretion through a magnetized 
disk during a T-Taur i star (TTS) phase. Most recently, 
iMatt fe Pudritj (|2008f ) report on 2-D MHD simulations of 
spindown by an aligned dipole field in a coronal wind during 
the TTS phase. With some translation for differences in no- 
tation and parameter definition, many of their results seem 
quite analogous to those reported here for hot-star winds. 
In particular, their eqn. (14) for the Alfven radius is quite 
similar to the dipole scaling in eqn. (I19|l here, just replacing 
Woo with the (numerically comparable) surface escape speed 
Vesc, with their numerical best-fit exponent m = 0.223 very 
similar to our 1/4 = 0.25. 



5 SUMMARY AND FUTURE OUTLOOK 

This paper analyzes the nature of angular momentum loss by 
radiatively driven winds from hot-stars with a dipole mag- 
netic field aligned to the stellar rotation axis. It applies our 
previous MHD simulation study from paper 2, which con- 
sisted of a 2-parameter series of models, dependent on the 
critical rotation ratio W , and on the magnetic confinement 
parameter r/ t . Key results can be summarized as follows: 



• As shown in figure[6] for both slow (W = 1/4) and mod- 
erate rotation (W = 1/2), the time-averaged, total angular 
momentum loss rate, J, follows closely a dipole-modified 
version of the Weber-Davis scaling of eqn. |T}. 

• Specifically, in the dipole case, the Alfven radius has 
a weaker dependence on magnetic confinement parameter, 
with a strong confinement scaling as Ra ~ riV A vs. as Ra ~ 
y/rjZ for the monopole [cf. eqns. (fT9)l and ©]. 

• This leads to a dipole scaling for J, given by eqn. (|20p . 
that is weaker than for a monopole, by a factor l/^/rj7 in the 
strong-confinement limit. It also implies a correspondingly 
longer spindown time. 

• As in the WD67 monopole case, the total angular mo- 
mentum loss is generally dominated by the magnetic compo- 
nent. However, in the strongest confinement models, there 
is trend toward increasing contribution by the gas, appar- 
ently a consequence of the eventual breakout of equatorially 
trapped material. 

• The dipole nature of the field gives the angular mo- 
mentum transport a complex variation in radius, latitude, 
and time, consisting of long intervals of gradual build-up, 
punctuated by episodic breakouts of material trapped in an 
equatorial disk. 

• However, the gradual buildup and storage of angular 
momentum in the circumstellar field and gas implies that 
the stellar spindown rate is likewise gradual, with no sud- 
den jumps during intervals of breakout. As such, "magne- 
tospheric emptying" does not seem like a likely explanation 
for sudden jumps in rotation rate claimed in some magnetic 
hot-stars. 

• The "dead zone" of closed loops surrounding the equa- 
tor does inhibit the equatorial loss of angular momentum 
from the near the stellar surface, but the net effect is merely 
to shunt the angular momentum flux to a magnetic compo- 
nent at mid-latitudes. The upshot is that the overall dipole 
scaling for total J effectively already accounts for any dead- 
zone reduction. 



The 2-D, time-dependent models here for rotation- 
aligned dipole thus provide a much richer physical picture 
for angular momentum loss than inferred in the 1-D, steady 
models for the WD67 idealization of a monopole field. How- 
ever, proper interpretation for the broad population of hot- 
stars inferred to have tilted dipoles or even multipole fields 
will require even more challenging 3-D MHD simulation 
models that explicitly allow for variations in azimuth as well 
as latitude and radius. Developing such fully 3-D MHD sim- 
ulations for these cases is thus a major focus of our planned 
future research. 
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APPENDIX A: ANGULAR MOMENTUM FLUX 
FROM MAGNETIC POLE FOR OBLIQUE 
ROTATOR 

To gain insight into the relative effectiveness of angular mo- 
mentum loss for the case of non-aligned dipole, let us ex- 
tend the WD67 monopole analysis to the 1-D flow directly 
along the magnetic pole of an oblique dipole, i.e., with axis 
now perpendicular to the rotation. This flow tube has a ra- 
dial orientation, but its areal expansion is some factor h(r) 
faster than r 2 , so that both radial mass flux and radial field 
st rength now vary as pv r ~ B r ~ 1/hr 2 . Following eqn. (20) 
of lOwocki fc ud-Doulal {2004), we can approximate this ex- 
pansion factor by 



, , . r R* + Ra 
Kr) = -- 



(Al) 



R* r + Ra 

where we take their "confinement radius" R c = Ra- 

The analysis then proceeds analogously to the WD 
monopole case summarized in section 2.3, leading again to 
an equation of the form (|14p . except that now the quantities 
rh = An pvr 2 , and thus j eq = J eq /rh, are no longer constant 
in radius. Instead, in terms of the constant mass loss rate 
M, we have rh(r) = M /h(r). Continuity at the Alfven criti- 
cal point again requires j eq = Q,R 2 A , but this now implies an 
angular momentum loss that scales as 



tlA 



(A2) 



where h A = h{R A ) = (l+R A /R*)/2. 

Since pv r ~ B r ~ 1/hr 2 along this dipole-axis flow, the 
ratio of magnetic to flow energy follows the scaling 

B 2 Ari.Rl 



r](r) 



(A3) 



4npv 2 h(r)r 2 (l- R^/r)! 3 ' 

where the factor 4 accounts for the fact that rj^ is defined 
in terms of B 2 q — Bp/ 4. Ignoring for simplicity the radial 
variation of velocity (effectively using /3 = 0), the Alfven 
condition r/(RA) = 1 thus implies (Ra/R*) 2 = &r)*/h,A, and 
so 

16r?* 



J e „ = MQRl^t- = MQR'i- 

hi 



(1 + Ra/R*Y 

In the strong confinement limit 77, 2> 1, we find Ra/R 
2r/l^ 3 , which gives 



(A4) 



J eq w M QRt 4ri, 



1/3 



rj* > 1- 



(A5) 



Comparing this to the monopole scaling J ~ r/* and the 
aligned-dipole scaling J ~ ^/rj*, we see that this oblique 

dipole-axis has an angular momentum loss rate that is 

1/3 

weaker than either of them, scaling only as J ~ r\J . 



